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Abstract 

We report the first observation of a near-threshold enhancement in the D^D^n^ system from 
B — ^ D^D^TT^K decays using a 414 fb~^ data sample collected at the T(45) resonance. The 
enhancement peaks at a mass M = 3875.4 it 0.7^20 MeV/c^ and the branching fraction for events 
in the peak is B{B D^W7r^K)={1.27 ± 0.31 toU) x lO""^. The data were collected with the 
Belle detector at the KEKB energy-asymmetric e+e" collider. 



PACS numbers: 13.25.Hw, 13.39.Mk, 14.40.Gx 



Belle recently discovered [l| a new state, the X(3872), with a mass of (3872.0 ± 0.6 ± 0.5) 
MeV/c^ and width less than 2.3 MeV/c^ in the J/t/^tt+tt^ system from — > J /ip'K^'K'K'^ 
decays. Other experiments confirmed the existence of the X(3872) 0], but all published 
results are in the J/ipTi^Ti~ mode only. Although the initial expectations were that the 
X(3872) was one of the unobserved charmonium states, subsequent experimental observa- 
tions disfavor this hypothesis [3]. 

Since the properties of the X(3872) are not consistent with a charmonium assignment, 
there have been speculations that it is some type of exotic state, for example, a qqg hybrid 
0. The X(3872) mass is within errors of the _D°D*° threshold (3871.2±0.9 MeV/c^), which 
triggered speculation that it might be a D^D*^ bound state (deuson) j^. If the X(3872) 
is a loosely bound S-wave molecule composed of D^D*^ charm mesons, it is expected that 
there will be an enhancement in the near-threshold D^D*^ invariant mass distribution 0,0]. 
The Belle collaboration found no evidence for D^D^, D^D~ and D^D^n^ decays of X(3872) 
with a smaller sample of BB events . The distributions of daughter particle momenta for 
D^D^TT^ decays of a D^D*^ molecule are expected to be different from those of an incoherent 
sum of the decays of free D*^ and 

In this letter, the D^Wti^ system is studied in B+ D°Wn'^K+ and 5° D^Wn^Kl 
decays. Inclusion of charge conjugate modes is implied throughout this paper. 

These results are based on a 414 fb~^ data sample corresponding to 447 million BB pair 
events collected with the Belle detector at the energy-asymmetric e^e~ collider KEKB 
The fractions of neutral and charged B mesons produced in the T(4S') peak are assumed 
to be equal. 

The Belle detector is a general purpose magnetic spectrometer with a 1.5 T magnetic field 
provided by a superconducting solenoid. Momenta of charged particles are measured using a 
silicon vertex detector and a 50-layer central drift chamber (CDC). Photons are detected in an 
electromagnetic calorimeter (ECL) consisting of 8736 CsI(Tl) crystals. Particle identification 
likelihoods, Ct^/k-, are derived from the information provided by an array of 128 time-of- 
fiight counters, an array of 1188 silica aerogel Cerenkov threshold counters and dE/dx- 
measurements in the CDC. 

Kaon candidates are selected from well-measured tracks by using a requirement on the 
likelihood ratio, Ck/{Ck + -^tt), which has an average kaon identification efficiency ~ 97% 
with a pion misidentification rate of ~ 18%. Similarly, charged pions are selected with an 
efficiency of ~ 98% and kaon misidentification rate of ~ 12%. All tracks compatible with 
the electron hypothesis (~ 0.2% misidentification rates from pion/kaon) are eliminated. 

Neutral pions (vr^) are reconstructed from pairs of isolated ECL clusters (photons) with 
invariant mass in the range 119 MeV/c^ < M^^ < 150 MeV/c^ (~ ±3(j). The energy 
of each photon is required to be greater than 30 MeV in the barrel region, defined as 
32° < 9j < 128°, and greater than 50 MeV in the endcap regions, defined as 17° < 6^^ < 32° 
or 128° < 6-y < 150°, where 6^ denotes the polar angle of the photon. A mass-constrained 
fit is applied to obtain the four-momentum of a ir^ candidate. 

Neutral kaons (Kg) are reconstructed via the Kg —>■ vr+vr" decay mode. There 
is no particle identification requirement for daughter pions and the requirement 
I M-n-TT — Mj^o \ < 11 MeV/c^ (~ 3.5(7, where a is the n~^7r~ invariant mass resolution) is ap- 
plied. Selection criteria to reduce random combinations of two tracks are described in detail 
elsewhere A mass- vertex constrained fit is performed to the Kg candidate to improve 



the resolution on its momentum measurement. 



Beam-gas events are rejected using the requirements |P^| < 2GeV/c and 1.0 < 
i?vis/-Ebeam < 2.5, where Pz, -Ebeam and i^vis are the longitudinal momentum sum, beam 
energy and total visible energy, respectively, in the center of mass (CM) frame. Continuum 
events (e"'"e~ — >• qq, where q = u,d, s, c) are suppressed by requirements on the ratio of the 
second to the zeroth Fox- Wolfram moments R2 < 0.50. 



Candidate mesons are reconstructed from well- measured charged tracks in the K~7[~^, 
K^TT^7T~^7r~ , Kg7r~^7^~ and K^K~ decay modes. A ±3 sigma mass window is applied for 
selecting D^'s, where sigma is the decay- mode-dependent resolution of the reconstructed 
mass (typically ~ 4.5 MeV/c^). Mass and vertex constrained fits are applied to improve the 
meson momentum resolution. 

A D^D^ candidate pair is combined with a 7r° and a kaon to reconstruct a B candi- 
date. Continuum events are further suppressed with the criterion, | cos6'Thrust| < 0.9, where 
^Thrust is the angle between the thrust axis of the B candidates and the thrust axis of 
the remaining tracks and isolated ECL clusters. The beam-energy constrained mass, Mbc 



(-Ebeam)^ — (Si-Pj)^)) whcre Pi is the momentum of the ith daughter of the candi- 
date B in the CM frame is restricted between 5.273 GeV/c^ and 5.286 GeV/c^. A peak in 
the difference between the measured energy of the B candidate and the beam energy, AE 
(=^j Ei — -Ebeam) is a siguaturc of i?-meson signal events, where Ei is the CM energy of the 
ith daughter of the candidate B. 

The B — i> D^D^TT^K signal Monte Carlo (MC) sample is generated in two steps, B 
X(3872)if followed by X(3872) — D^D^n^ assuming a phase space distribution in both 
decay chains. The average number of D^D^ti^K entries per MC signal event is ~ 3.55, 
which are mainly due to multiple slow tt^'s. The characteristics of incorrectly reconstructed 
7r° candidates are identical to signal tt^'s. To reduce this multiplicity we use the selection 
criterion, M£)0^o < 2.013 GeV/c^ or iV%o^o < 2.013 GeV/c^. This requirement reduces 
background and the candidate multiplicity to 1.68 with almost no loss of signal efficiency. 
Possible bias due to this selection criterion is studied in the following event samples: (z) a 
large sample of generic BB and continuum MC events; (u) D^D-ti^K data; {in) D^D^n^K 
(same-flavor charm) data; {iv) the D°-mass sideband data (one Z}°-meson is reconstructed 
when the invariant mass of daughters is outside the D°-mass signal region); and (f) AE 
side-band data (60 MeV < \AE\ < 110 MeV). No peaking behavior is observed in the DDn^ 
mass distribution for any of the above-mentioned control samples, thereby confirming that 
there is no bias in the selection criteria on M^otjo^o. 

A unique D^D^tt^K candidate is chosen out of possible multiple candidates in a given 
event by taking the combination with the minimum value of ( ^) + ( ^) + 



where Ax and cr^ are the deviation of the measured quantity x from its nominal 

value and the uncertainty in its measurement, respectively. Multiple kaon entries are re- 
solved by choosing the candidate with the highest kaon identification probability for charged 
kaon and minimum {Mt^t^ — Mj^o\ for neutral kaon. There is a negligibly small number of 
events with charged and neutral kaon multiple entries. 

An unbinned extended maximum likelihood fit to the D'^D^n^ invariant mass, M^o;do^o, 



and AE distributions is used to obtain the signal yield. The fit includes three components: 
(i) a signal function, which is modelled by the sum of two Gaussian functions with the 
same mean value for AE and single Gaussian function for M^^otjo^o; {ii) a non-resonant 
B D^D*^K signal, where AE is also modelled by a double Gaussian function and M^otjo^o 
with a threshold function; and {Hi) the remaining backgrounds, which are modelled with a 
first-order polynomial for AE and another threshold function for Mj^o-^o^o, this threshold 
function is obtained from the AE sideband data oi B ^ D'^D^tc^K events. Shapes of the 
AE distributions for signal and non-resonant B D^D*^K background are fixed from 
the B D^D*^K data sample. The signal has a narrow Gaussian component with width 
a ~ 4.5 MeV and a wide Gaussian component with width 4.6 times larger that accounts for 
40% of the signal. 

Parameters of the M^otjo^o threshold functions are fixed from a large MC sample of 
B D^D*^K events for the non-resonant components and B D'^D^tt^K AE sideband 
data for remaining backgrounds. The normalization factor for the non-resonant component is 
fixed according to the branching fraction from 3| • The slope of the background polynomial. 



the parameters of Mj^ojjo^o threshold peak, and the normalization factors of signal and 
combinatorial background component are free parameters of the fit. 

Figure Ufa) shows the scatterplot of AE and M^oj^Ot^o in data. There is a cluster of events 
in the D^D^n'^ threshold region. The AE distributions for different M^otjo^o intervals are 
shown in plots (b)-(m), where a one-dimensional fit gives a signal of (23.4±5.6) events in 
the M^o7)0^o range from 3.870 GeV/c^ to 3.878 GeV/c^. The statistical significance of this 
signal, defined as ■>/— 2 In (£o/'Cmax), where £o(max) is the likelihood without (with) the 
signal contribution, is 6.4(T. A similar analysis that uses the Mbc distribution rather than 
the AE distribution to measure the signal also shows a clear peak for the same M^o-^o^o 
interval, with a consistent signal yield and similar statistical significance. 

To obtain the exact position of the near-threshold peak as well as its branching fraction, 
the two-dimensional distribution of AE and Q-value (=^^0750^0 — 2M£)0 — M^o) is fitted. 
Projections onto Q- value (for \AE\ < 25 MeV) and AE (for 6 MeV/c^ < Q- value < 14 
MeV/c^) are shown in Figure |21 along with the results of the fit. The fitted mean and width 
of the near- threshold peak are 11.21±0.68 MeV/c^ and 2.42±0.55 MeV/c^ respectively, 
in the Q-value. The signal yield is 24.1±6.1 and the significance, including the effects of 
systematic errors, is 6.4(T. Individual results for the charged and neutral B meson samples 
are given in Table E] together with the combined result. 

In terms of the invariant mass of D^D^tt'^ system, the peak position is M^oj^o^o = 3875.4 
±0.7 MeV, where the error is statistical only. 



TABLE I: Signal yield and branching fraction for B —>■ D^D^tt^K with near-threshold peak near 
3.8754 GeV/c2. 

Signal eBxW^ Nobs sig, a B xW^ 

OTiO^O TV n io_Ln in O/I 1 j-« 1 a a 1 ov^ n 01 +0.22 



-0.39 



B D^D'^-K^K 2.12±0.10 24.1±6.1 6.4 1.27± 0.31 
B+ DW^'k^K+ 3.62±0.14 17.4±5.2 5.0 1.07± 0.3lj:[!:;^^ 
50 ^ d^Wtt^K^ 0.84±0.04 6.5±2.6 4.6 1.73± 0.70toil 



The MC-determined signal efficiency for the near-threshold peak is (1.87±0.05)%. The 



3.95 F 



O 

o 

O 

Q 



O 



o 



o 



I 
o 



|:|,ga=:Bia:g§B:=|:=,^: :bB::s 



3.9 



ig§sTiisi^ill 

==ii:„e9-=i-- 



(a) 



-0.1 



AE (GeVl 




20 



3.& 2<M 



i 



:3.870 




28 



(e) 3.886<M<3.894 




30 



(h) 3.910<M<3.918 




(k) 3.934<M<3.942 




(1) 3.942<M<3.950 



0.1 

16^ 



0.2 




24 



(d) 3.878<M<3.886 



I 





32 



(g) 3.902<M<3.910 




ill 





46 



(]) 3.926<M<3.934 




(m) 3.950<M<3.958 



T 




0.2 
A E (GeV) 



0.2 
A E (GeV) 



0.2 
A E (GeV) 



FIG. 1: (a) Scatter plot of AE and Mj^q-^q^q distribution in the data sample, (b)-(m) AE distri- 
butions in 8 MeV/c^ D^D^-k^ mass bins for a possible B —>■ D^D^ir^K signal in data. 



contribution from non-resonant B — > D^D*^K events in the near-threshold Mj-)0-j^o^o region, 
calculated from a large sample (~ 60 times real data) is found to be 1.6 ± 0.2 events. 

Because of the Hmited available phase space, the Moo^o/A%o^o distributions near the 
threshold show some clustering around D*" mass. Thus it is not possible to separate the 
contributions of D^D*^ and D^D^tt^ to the peak. 

With a large sample of MC-simulated D^D^n^ events generated in a narrow, near- 
threshold mass peak, it is found that the reconstructed peak in the M^o-po^o distribution has 
a high- mass tail that is caused by poorly reconstructed tt^'s. We are unable to distinguish 




FIG. 2: Projection of the unbinned fit (a) Q-value {=Mj-)0-j^Ot^o — 2Mj^o — M^a) in /S.E signal region 
(|A£;| < 25 MeV) and (b) l^E in the signal region of Q-value (6 MeV/c^ < Q-value < 14 MeV/c^). 
The dots are data points, the hatched histogram corresponds to combinatorial background; the 
dashed line indicates the total background and the solid line is from the combined fitting function. 



this high-mass tail component in the data or in samples of MC signal plus background with 
sizes similar to the data. This tail, if it exists, would produce a positive bias on the peak 
mass measurement. We account for this possibility with asymmetric systematic errors on 
the peak mass and efficiency of l^'y MeV/c^ and 1^|'8%, respectively. Including all sys- 
tematics, the observed M^otjo^o peak position is 3875.4 it 0.7;!^?;7 ± 0.9 MeV/c^, where the 
second error is mainly due to the calibration uncertainty of the vr" energy and the effects of 
a possible high-mass tail as discussed above. The third error is due to the uncertainty in 
the world-average mass [l^ . 

The systematic uncertainty on the B D^D^'k^K branching fraction for the near- 
threshold peak is obtained from the quadratic sum of the uncertainties due to (a) limited 
MC statistics (1.3%), (b) subtraction of the B D^W^K contribution (5.1%), (c) num- 
ber of BB events {Ng^) (1.3%), (d) PDG branching fraction of and K^g (5.0%), (e) 
track finding efficiencies (9.5%), (f) i^/vr identification uncertainties (7.0%), estimated us- 
ing D*~ —>■ D°(— > K^TT~)7r~ events, (g) 7r° detection efficiency (7.0%), estimated from a 
comparison of K^tt^tt^ yields in data and MC, (h) Kg selection efficiency, estimated 

from a comparison of /C^tt+tt^ yields in data and MC (2.1%), (i) the ratio of D°-mass 

window in data and MC (2.0%), (j) signal efficiency, calculated from the difference in AE 
and Mj^ojjo^o fits (5.8%), (k) efficiency due to poorly reconstructed 7r° (125.9%) 0) ^lass 
value (1.0%). The total uncertainty is estimated to be 

In summary, a near-threshold D^D^n^ invariant mass enhancement is observed at 
3875.4±0.7l'i':7±0.9 MeV/c^ in S ^ D^Wtt^K decays. The significance of this enhancement 
is 6.4(T. 

The observed D^D^tt^ mass is 2.0o" higher than the world-average value of the X(3872) 
mass of 3871.2 ±0.5 MeV/c^ while the branching fraction of this threshold peak is 9.4t^ 
times larger than B{B+ X{3872)K+) x S(X(3872) J/ipn+n-) Q. Reference |l3 
ruled out all possible quantum states of X(3872) except J^^ = and 2+"*" while CDF 



finds that possible quantum number assignments are l^"*" and 2 + [16]. If this near-threshold 
enhancement is due to the X(3872), the J'^'" = l^"*" quantum number assignment for the 
X(3872) is favored, because the near-threshold decay X(3872) D^'D*^ / D^W'k^ is ex- 
pected to be strongly suppressed for J = 2. 
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